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We consider a two-dimensional spin system in a honeycomb lattice configuration that exhibits 
anyonic and fermionic excitations [Kitaev, cond-mat/0506438 . The exact spectrum that corre- 
sponds to the translationally invariant case of a vortex-lattice is derived from which the energy of 
a single pair of vortices can be estimated. The anyonic properties of the vortices are demonstrated 
and their generation and transportation manipulations are explicitly given. A simple interference 
experiment with six spins is proposed that can reveal the anyonic statistics of this model. 

I. INTRODUCTION 

Lately, a remarkable eflFort has been made towards the understanding of topological systems (see e.g. 0, 0, ll| 
j^"*' and references therein) . This ranges from the realization of systems that exhibit topological behavior 0, 0, |0| to the 
systematic identification and characterization of their topological properties 0, S lol [lo|. [TH IT^ . These efforts 
have been greatly boosted by an interest in performing error-free quantum computation iTHlsf. The idea is to take 
\^ [ advantage of the non-trivial statistical properties of anyonic particles, that deviate from the bosonic or fermionic 
behavior, in order to encode and process information. In these systems, quantum computation is protected from any 
local errors that do not destroy the nature of the anyons. 

This article considers in detail a specific two dimensional lattice model presented by Kitaev It comprises of 
spin- 1/2 particles in a honeycomb configuration that are subject to nearest neighbor interactions. The corresponding 
quadratic Hamiltonian has been solved analytically in the vortex-free sector and a perturbation treatment was given 
for its vortex sector Here we present a non-perturbative study of the vortex-lattice sector. The latter consists of 
. vortices placed at each hexagonal plaquette of the lattice. From these analytic results one can estimate the spectrum of 
^3 ' two individual vortices in the limit that the interactions between vortices are negligible. Our findings are supported by 
[ numerical diagonalization of a small lattice system with 16 spins. A detailed presentation is given of the manipulation 
I procedures required to generate or transport the anyons. The considered system can be employed to perform robust 
^ i' quantum computation with abelian anyons [3 | . 

I ' , , This article is organized as follows. In Sectional we present the model and the symmetry properties that enable 
P5 ' its analytic treatment. In Section [llll the separation of the spectrum into different vortex sectors is presented and the 
^ , properties of anyonic vortices are elaborated. The particular cases of the vortex-free and vortex-lattice sectors are 
considered and their spectrum is explicitly derived. In the conclusions, the properties of the spectrum are analyzed 
and an exact numerical treatment is presented that supports our findings. Finally, a minimal lattice cell with six spins 
is presented that supports anyonic excitations. A simple interference setup is proposed that can reveal the anyonic 
statistics experimentally. 
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II. PRESENTATION OF THE MODEL 



Consider a two dimensional system with spin-1/2 particles located at the vertices of a honeycomb lattice The 
spins are assumed to interact with each other via the following Hamiltonian 

x-links y-links z-liiiks 

where "x-links" , "y-links" and "z-links" are depicted in Figure ^ Here we shall adopt periodic boundary conditions 
although open boundary conditions can be similarly treated. Consider an individual hexagonal plaquette and the 
corresponding operator 



(2) 



where af is a Pauli operator at vertex i of the hexagon. This operator commutes with the Hamiltonian, H , as well 
as with the corresponding operators, Wp' , of all the other plaquettes. Thus, its eigenvalues, Wp = ±1, are conserved 
quantities. Hence we can split the Hilbert space in sectors, characterized by a fixed configuration of Wp's and 
solve for the eigenvalues of the Hamiltonian in each sector. A plaquette with Wp = —1 corresponds to a vortex with 
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anyonic statistics as we shall see in the following. Here the exact diagonalization of the Hamiltonian will be presented 
for the limiting cases, where Wp — 1 or Wp — —1, for all plaquettes, p. They correspond to the vortex-free and a 
regular vortex-lattice configuration, the latter having a vortex at each hexagonal plaquette. 

To diagonalize Hamiltonian we rewrite the spin operators in terms of Majorana fermions. The latter are defined 
as the "real" and "imaginary" parts of usual fermionic operators, and a^, in the following way 

I flfe — a! 

C2k-i = cik + al, C2fe = -. . (3) 

The Hcrmitian operators, q, satisfy the following relations, cj = Ci and {ci,Cj} = 26ij. Here we introduce the 
fermionic picture by representing each spin with two fermionic modes, oi and 02 in such a way that the no-fcrmion 
state represents a spin up and the two-fermion state represents a spin down. There are four Majorana operators that 
correspond to each spin given by 

t '^i ~ 'A 1 X , t 7iy a2 — at 

Co = fli + a = c, ci = — b , C2 — a2 + a!, = 0'', C3 — — b . (4) 

i " i 

It is possible to project the full space of states of the two fermions to the physical space C of the spin states by 
employing the projection operator D — b^b^-'b'c, where 

I *> e £ ^ D I = I . (5) 

Within the subspace C the following identifications hold 

a" — ib°'c, {a = x,y,z), (6) 



where a" is a specific representation of the Pauli algebra. 

To employ this representation for Hamiltonian ^ we introduce four Majorana fermions 5| and Cj for each site 

j. As [a", Dk] = the diagonalization of Hamiltonian ^ is re-expressed as diagonalization of H{{a"}) = H{{ib"cj}) 
accompanied by the constraint Dj | = | ^'). With the above transformation we have that <^'j<^'k = ^iujkCjCk with 
Ujk = ib'^b'^, finally obtaining 

* V i r r, where 4 , - / ^-^a^jfe, if j and k are connected by a link, . . 

^ - ll^^jkCjCk, where A,k - < otherwise. 

Observe that [H,Ujk] — 0, ujj. = iijk and Ujf. = 1, hence, one can restrict to an eigenspace of states given by a 
certain configuration of eigenvalues (ujk — ±1) of the Ujk operator. This reduces the diagonalization of Hamiltonian 
Q to the diagonalization of the corresponding Hamiltonian, where Ajk is substituted by a certain configuration Ajk 
obtained by replacing Ujk with its chosen eigenvalue. Thus, instead of considering the decomposition of the Hilbert 
space with respect to eigenvalues of lip one can consider the decomposition, £„, with respect to the eigenstates of 
Uij. Restricted in the physical subspace C the following relation holds 

Wp= Y\_ (8) 

{j,k)£ houndary{p) 

where the (j, k) links are ordered in a clockwise fashion around the plaquette boundary. However, the operators Ujk 
and Dj do not commute so the resulting eigenstates, | of the Hamiltonian do not satisfy the constraint This 
is remedied by the symmetrization 

I = n (^^) I *«) ■ (9) 

Thus, the eigenstates of Hamiltonian Q have to be symmetrized with the above procedure in order to obtain the 
physical states that correspond to the vortex configurations of Hamiltonian Obviously this procedure does not 
affect the eigenvalues of the Hamiltonian. 
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III. SPECTRUM OF HAMILTONIAN 



A. Sectors of Hilbert space 

As we have seen, it is possible to reduce the full Hilbert space of Hamiltonian |^ to independent subspaces, C^, that 
correspond to different eigenvalue configurations of Wp. The subspace with Wp — +1 for all plaquettes p corresponds 
to the vortex-free configuration. Diagonalizing Hamiltonian ((T)) in this subspace will just give a fermionic spectrum 
with no anyons (vortices) present. Changing the sign of one Ujk results in two adjacent plaquettes having Wp = —1. 
This configuration corresponds to two vortices placed at these plaquettes. In the absence of external fields these 
anyonic excitations are static in the sense that they remain there indefinitely. One can perform this change of sign in 
Ujk by simple Pauli rotations of the original spins. Indeed, a rotation on a certain spin will cause the generation 
of two vortices in the adjacent plaquettes, as shown in Figure ^ Hamiltonian Q can now be diagonalized in this 
subspace of states which results into a fermionic spectrum. The ground state of each sector corresponds to the pure 
anyonic configuration without fermions. 





z z 


z 














FIG. 1: The honeycomb lattice with the three different types of links. At each type corresponds a different spin-spin interaction. 
Two vortex configuration is induced by a rotation at a certain spin. This rotation does not commute with the interaction 
terms on the x- and j/-links of Hamiltonian a property represented by the dotted connecting string. 



To study the vortex excitations as well as their transport properties systematically we shall adopt the operators 
that correspond to the toric code limit of the model 2, 14] obtained, e.g. when Jz ^ Jx,Jy- In this limit there are 
three vortex excitations above the ground state, | gs), given by 

I Z) = (7^ 1 1 5s) « 1 ® cr^ I <?s) , I r) = (j^ ® I gs) « (g) cr^ I gs) , \X)= ia"" (g) | gs) . (10) 

The two operators of the tensor product are acting on the two spins of a certain z-link. The | Z) and | Y) excitations 
behave as bosons with themselves while the | X) excitation is fermionic. The latter can be produced by fusing the 
I Z) and I Y) particles and its energy is the sum of the energies of the two constituents. Appropriate applications of 
the Pauli operators result in the transport of the anyonic subspaces around the lattice. All of the excitations behave 
as anyons with each other with a statistical angle 6 — tt/2 that results from the anticommutation relations of the 
Pauli operators. Similarly, the fusion properties of the vortices follow from the properties of the Pauli operators, 
XxX = YxY = ZxZ^l, XxY = Z and permutations, where 1 is the no- vortex configuration. For general 
values of the couplings Ja these excitations may be accompanied by non-vortex fermions that reside at each sector. 
This, of course, does not affect the anyonic properties of the vortices. Hence, it is possible to move between different 
subspaces, Cw, by applying Pauli rotations on the original spins. 

The generation and transport properties of the different sectors do not necessarily correspond to the properties of 
the anyonic excitations without the fermions. In general, a simple rotation moves the vortex-free configuration to 
the one with two vortices, but fermions may also be generated. To guarantee that no unwanted fermions are created 
one may consider realizing locally a regime for which one coupling is much larger than the other two, e.g. Jz 3> Jx, Jy 
Then the fermionic energy gap is much larger than the anyonic one 'l4'|, while we can still find Pauli rotations that 
correspond to the exchange of energy of the order of the anyonic gap. Indeed, a rotation of spin i of the initial 
Hamiltonian results in a two vortex configuration 

i/2-v = <i?^f =H + 2 J,<aJ + 2JyC7^al (11) 

where (i,j) denotes an a;-link and {i,k) denotes a y-link, and, hence it causes an exchange of energy that scales 
with Jx and Jy. As we shall see in the following, the fermionic excitations are of the order of Jz- Hence, they are 
decoupled from the rest of the system and are negligible in manipulations at this energy scale. It is worth noting that 
diagonalizing Hamiltonian (|ll|l does not give a rotated | Z) state (which would correspond to the lowest energy 
vortex-free configuration) as that state would correspond to a different sector. 
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B. The vortex- free sector 



Due to the symmetry of the Hamiltonian we can restrict ourselves to the case of positive couphngs J a ,2] . The 
simplest configuration of the eigenvalues of Ujk is given by Ujk = 1 for all links (j, k) of every hexagonal plaquette 
taken in a clockwise fashion. This configuration is also known as the vortex-free configuration as it gives Wp = \ for 
all plaquettes p, and can easily be solved by a Fourier transformation. 
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(a) 



(b) 



(c) 



FIG. 2: (a) The effective square lattice of the original honeycomb model. Each cell (filled circle) incorporates two vertices of 
the original lattice and a z-link, while the depicted x and y-links have unit length, (b) The vortex-free case where Ujk = -1-1 
for {j,k) taken clockwise around each plaquette. (c) The vortex-lattice where Ujk = for all (j, fc) oriented clockwise except 
for the ones that correspond to every other z-link in the x direction. 



It is convenient to rearrange the lattice in the following way. Let a unit cell include the (j, k) site that corresponds to 
the z-link. Each cell can be represented by an index s while the position of the vertex within the cell (i.e. distinguishing 
between j and k) is represented by an index A (see Figure|2Ia)). Thus, the Hamiltonian takes the form 



As Hamiltonian l)12|l is translationally invariant we employ Fourier transforms to diagonalize it. Let us define 



(12) 



(13) 



where Px,y € tt] and the two dimensional vector, t = {t^, ty), indicates the position of the vertices on the effective 
square lattice. The Hamiltonian becomes 



where c^(p) = cx{-p), cx{p)cUq) + d{q)cx{p) = Spq6xf, and 



iA{p) = 



if{p) 

-irip) 



, f{p) = 2(J,e^f== + Jye'Py + J,). 



For = arg(/(p)) one can introduce the fermionic operators 

b{p) = (ci(p) + ie''^C2{p))/V2, b\p) = {ci{-p) ~ ie-''^C2{-p))/V2. 
In terms of b the Hamiltonian takes the diagonal form 



I d^\f{p)\ 



b\p)b{p) 



The ground state is given by 



(14) 



(15) 



(16) 



(17) 



= n ^^p) 1 0) = n (^i w + »e-'^4b)) 1 0) 



(18) 
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where 1 0) is the vacuum state of the Majorana operators, c\{p). The overall ground state energy is Egs = 
— i / cPp \f{p)\- The first excited state and its energy gap above the ground state are given by 

|epo> =&^(Po)|0), Ai?p„ =min|/(po)|. (19) 

Po 

It is clear that if \Jx\, \Jy\, \Jz\ satisfy 

\J.\<\Jy\ + \Jzl \J.\ < \Jx\ + \Jyl \Jy\<\JA + \Jxl (20) 

then there exists a po such that AEp^ — 0, i.e. the gap vanishes (see Figure IHl^a)). The excitation spectrum of the 
operators 6(po) corresponds to the fermionic spectrum of this sector (wp = 1 for all plaquettes, p). Assuming that 
the coupling configuration is not at the borders of the gapless region one can easily verify that in the gapless region 
there are two distinctive momenta for which the energy gap becomes zero. At these two Fermi points the energy has 
conical singularities. 




(a) gapless (b) 



FIG. 3: (a) The gapped and gapless phases for the vortex-free case, (b) The gapped and gapless phases for the vortex-lattice 
case. For both plots Jx + Jy + Jz ~ 1 where we assume Jx,Jy,Jz > 0. For the vortex-free case (a) the inner boundaries 
are given by the conditions Jx + Jy ~ Jz, Jy + Jz = Jx and Jz + Jx = Jy, while for the vortex-lattice (b) the conditions are 



C. The vortex-lattice sector 



The anyonic spectrum can be obtained by studying the configuration with Wp = —1 for all p. To implement 
it we consider the particular case where ujk has alternating signs at the z-links along the x direction while it is 
homogeneous in the y direction (see Figure I^Jc)). This is sufficient to create Wp = —1 for all plaquettes p giving rise 
to a vortex-lattice where a vortex is placed at each hexagonal plaquette. 

As the sign of the z-links coupling is now alternating compared to the previous vortex-free case, the Hamiltonian 
in the momentum representation takes the form 



H^'^Y^J d^d'p'A^,{p,p')cx{p)Cf,ip'), 



(21) 



where Ai2{p,p') = S^p + p')M{p) + 5{p + p' + Ti)N{p), M{p) = 2 J^^e^P- + 2Jye'Py , N{p) = 2J^ and A21 = -Al^- 
Explicitly we have 



U = 



dpxdpy [M{p)c\{p)c2ip) + M*{p)ci{p)c\{p) + N{p) [c\{p)c2{p) + ci{p)c\{p)\ } , (22) 



where p = {px,Py) and p — {px ± TT,py). As the last term of ()22(l is repeated twice in the interval px G [— 7r,7r] the 
Hamiltonian can be rewritten as 



H 



7r/2 
-7r/2 



dpx j dpy i^M(p)c\{p)c2{p) + M(p)c\(p)c2{p) + Nip) [c\{p)c2{p) + c\{p)c2{p)\ } 



h.c. 



(23) 
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To diagonalize the Hamiltonian ^3 '^^ would like to transform (|23|l into the following form 
H = t: dpx I dpy 

^ J-7T 

The corresponding canonical transformation has to satisfy 



A{p)4{p)r,2{p) + B{p)r,l{p)f]2{p)\ + h.c. (24) 



[H,r,l{p)] = -'-A*{p)4{p), [H,4{p)] = -'-B*{p)4{p), [H,r,lip)]^'-Aip)r^l{p), [H,4{p)] ^ '-Bip)rjl{p). 
One can verify that the appropriate transformation is given by 

ca(p) = cos0r]x{p) - sin6'e"'"^'^*'^77A(p), 

cxiP) = sin0e(-i)'*^r,A(p) +cos0?7a(p), (25) 

where </> = arg(M (p) + M*(p)) and 

tan2.= W^) + MW. (^6) 
|M(p)|2 - |M(p)|2 ^ ' 

The coefficients A{p) and B{p) of Hamiltonian (|24|l are given by 

|Af(p)|2 + |M(p)|2 + 2Af(p)Af(p)- 



A{p) 



M*{p) + M{p) 



|M(p)p + |Af (p)p + 2M{p)M{p) - VA 
^ A/(p) + Af*(p) ' 

where 

A = (|A/(p)|2 - |A/(p)|2)' + |A/(p) + M*{p))\'N\ (28) 
Finally, the Hamiltonian can be diagonalized by employing the following fermionic operators 

b{p) = (,71 (p) + ze'^,72(p))/V2, b{p) = {m{P) + ie'^'r72(p))/V2, 

fet(p) ^ (y^t^p) „ je-*'^^t(p))/^^ ;,t(p) ^ (,,t(p) „ ,e-^'^',,t(p))/y2. (29) 
where (p = arg(j4(p)) and ip' — arg(i?(p)), eventually giving 



/7r/2 />7r 
rfPa; / dpy 
-7r/2 J-TT 

The ground state with its energy is given by 



\A\b\p)h{p) + \B\b\p)b{p) 



\A\ + \B\ 



/2 . r . \A\ + \B\ 



(30) 



|gs)= II b{q)b{q)\Q), i^gV--/ rfP. / ^P, ' ' „ ' ' , (31) 
— 7r/2<g<7r/2 ' 

while some of the first excited states with their energy gap are given by 

I epo> = b^iPo) I gs> , AEp„ = min |^(po)|, 

Po 

I epo) = b\po) I gs) , AEp, = min \B{po)\. (32) 

Pa 

It is possible to derive the values of the couplings for which the vortex-lattice configuration becomes gapless. Note 
first that |v4(p)| cannot be zero, but |-B(p)| can. Indeed, there exists p such that |-B(p)| = if and only if all of the 
conditions below are satisfied 

which is given schematically in Figure Otb) . When H33|) is satisfied and for coupling configurations away from the 
borders of the gapless region one can verify that there are two distinctive momenta for which the energy gap becomes 
zero. At these two Fermi points the energy has conical singularities. 
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IV. CONCLUSIONS 



As we have seen, vortices with respect to the Ujk gauge field have anyonic statistics and are generated in the system 
by applying proper spin rotations. They are connected either with other vortices, or, when the system is open, with 
its boundary by string operators. The vortex excitations are static, but they can be transported by applying external 
fields. Each vortex configuration is accompanied by a fermionic spectrum. We evaluated this spectrum for the limiting 
cases of the vortex-free and the vortex-lattice configurations. Their ground states are fermionic vacua (pure anyons) 
from which one can obtain the physical wavefunctions by performing symmetrization Q. 

It is of interest to study the behavior of the energy of these configurations as functions of the couplings J a ■ When 
Jx = Jy = then Eg~^ = E'^l^ = —J^N where A'' is the total number of plaquettes for both the vortex- free and the 
vortex-lattice configurations. It is possible to check that for all other values of the couplings Jq we have Elf < E^f. 
Hence, the energy minimum is achieved by the vortex-free configuration, a fact that also follows from a theorem by 
Lieb 'l6l and has been verified numerically by Kitaev P|. It is easily observed that the two energies become equal 
only when one of the couplings Ja becomes zero. 

One can compare the fermionic excitation energy of the vortex-free configuration with the energy gap of a pair 
of vortices. To derive the energy of this pair from the energy of the vortex-lattice we assume that the vortices are 
non-interactive. This assumption is valid in the gapped regime where any interaction mediated by the fermions is 
exponentially suppressed |l7| . This is also supported by exact numerical diagonalization as we see in the following. 
In Figure ^a) the energy gap is plotted between the ground state and the anyonic or the fermionic excitations. One 
can easily see that there is a transition in the character of the first excitation from anyonic to fermionic. 



AE 
1.4-, 




W M 

FIG. 4: (a) Analytic results for the energy of a pair of anyons and of a fermionic excitation. The anyonic and fermionic energy 
gaps are plotted against Jx and Jy where we have set Jz — 1. (b) Numerical evaluation of the energy gap of the first excited 
state above the ground state. This gap is in agreement with the one given as a minimum of the vortex and fermionic gaps as 
seen in (a). 



To verify our findings we simulate the original Hamiltonian with 16 spins in a periodic honeycomb lattice. In Figure 
0fb) we plot the energy gap between the ground state and the first excited one against Jx and Jy, while keeping 
Jz = 1. There we see that the excitation becomes gapless for values of the couplings which agree with the behavior of 
the vortex and fermionic excitations. Indeed, the energy gap obtained numerically is approximately nested below the 
analytically obtained energy gaps with a maximum value that agrees well with the value obtained from Figure 2Ja). 

We conclude by presenting a simple model with six spins that can reveal the anyonic statistics of the excitations 
experimentally. It consists of one hexagonal plaquette with only the spin interactions that lay on the hexagon (see 
Figure [SJ. An anyonic excitation in the plaquette can be produced by a dg rotation. Its paired anyon is outside our 
system thus it does not increase the total energy of the system. In order to obtain the anyonic properties one would 
like to circulate another anyon around it. Even if it is not possible to have two anyons present in this simple model, 
one can consider anyons outside the system that can circulate the existing anyon by appropriately rotating the six 
spins. Indeed, performing the operation S = Ciclcrlcrlcrgcrl corresponds to a looping trajectory around the anyonic 
excitation in the plaquette and, thus, it should generate an overall tt phase factor. Indeed, for the ground state, | gs), 
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and the vortex state, | v) 




\v) 



(34) 



The property S \ gs) — \ gs) is due to the invariance of the ground state under closed loop operations. It is worth 
noting that S is also the Wp operator 101 for the single plaquette of our system. 



FIG. 5: A hexagon that can support anyonic statistics. A vortex is generated, e.g. by a rotation of spin 6. Pauli rotations 
around the plaquette result in the generation of the looping trajectory. This procedure results in a tt phase revealing the 
anyonic statistics of the vortex. 



This procedure can be used in an interferometric setup to observe the anyonic statistics experimentally. Indeed, if 
one performs the rotation R = {1 — ztT|)/\/2, then the superposed state (| gs) — i \ v))/\/2 is produced. The looping 
evolution and the inverse rotation will result into the | v) state if the statistics of the vortex excitation is anyonic, 
while a bosonic or fermionic statistics would give back the | gs) state. This interference effect produced with only one 
hexagonal plaquette can, in principle, be realized with present technology using NMR ^3 or Josephson Junctions 
Alternatively, one can consider a lattice of vortices where every other vortex is transported around its static 
neighbor. Such an experiment could be realized in an optical lattice setup with atoms |23| or polar molecules The 
readout of the final configuration could be performed by energy addressing using the spectrum derived in Section [llll 
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